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Abstract: The Chinese tree shrew (Tupaia belangeri chinensis) is a small experimental animal with a close affinity to primates. This 
species has long been proposed to be an alternative experimental animal to primates in biomedical research. Despite decades of study, 
there is no pure breed for this animal, and the overall genetic diversity of wild tree shrews remains largely unknown. In order to 
obtain a set of genetic markers for evaluating the genetic diversity of tree shrew wild populations and tracing the lineages in 
inbreeding populations, we developed 12 polymorphic microsatellite markers from the genomic DNA of the tree shrew. An analysis 
of a wild population of 117 individuals collected from the suburb of Kunming, China, showed that these loci exhibited a highly 
expected heterozygosity (0.616). These 12 microsatellites were sufficient for individual identification and parentage analysis. The 


microsatellite markers developed in this study will be of use in evaluating genetic diversity and lineage tracing for the tree shrew. 
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The Chinese tree shrew (Tupaia belangeri chinensis) 
is a squirrel-like animal. It belongs to the Family 
Tupaiidae of Scandentia. Tree shrews have the highest 
brain-to-body mass ratio of known mammals (Peng et al, 
1991). Because tree shrews share some characteristics 
with primates and insectivores, the exact phylogenetic 
position of the tree shrew has been debated (Arnason et 
al, 2002; Janecka et al, 2007; Peng et al, 1991; Xu et al, 
2012; Xu et al, 2013a). Analyses of the mitochondrial 
DNA (mtDNA) genome have shown that the tree shrew 
has a closer affinity with Lagomorpha (Arnason et al, 
2002; Xu et al, 2012), while nuclear gene sequences 
provide evidence for a close affinity to primates (Janecka 
et al, 2007; Killian et al, 2001; Lindblad-Toh et al, 2011). 

The tree shrew has long been proposed to be an 
alternative experimental animal to primates in 
biomedical research due to its characteristics such as 
small body size, short reproductive cycle and life span, 
and low-cost maintenance (Peng et al, 1991). This animal 
has been used to create animal models for infection with 
the hepatitis B and C viruses, the development of 
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hepatocellular carcinoma, myopia, psychosocial stress 
(Cao et al, 2003; McBrien & Norton, 1992; van Kampen 
et al, 2002; Yan et al, 1996a; Yan et al, 1996b) and 
learning and memory (Wang et al, 2011), although there 
are many problems awaiting further studies (Xu et al, 
2013b) . Despite an increasing interest in using the tree 
shrew to establish animal models for medical and 
biological research, most of the animals used were 
captured from the wild and/or domesticated for a few 
months. There is no pure inbreeding strain with a clear 
genetic background similar to that of a mouse or rat at 
the present time, which would largely influence the 
stability and repeatability of experiments. 

Tupaiidae contains four genera: Tupaia, Anathana, 
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Urogale, and Dendrogale (Helgen, 2005). Currently, a 
total of 15 species is recognized in Tupaia, which is 
broadly distributed across South and Southeast Asia, 
including southern China, India, Philippines, Java, 
Borneo, Sumatra and Bali (Olson et al, 2005; Peng et al, 
1991). Chinese tree shrews (Tupaia belangeri) are 
distributed across southwest China and are divided into 
six subspecies according to geographical distribution and 
morphological characteristics (Wang, 1987). 
Understanding the genetic diversity of the wild tree 
shrew population has critical importance for 
conservation and for a breeding program to establish 
inbred lines. Recently, we evaluated the mtDNA genetic 
diversity of the tree shrew inhabiting urban Kunming and 
observed a relatively high diversity (Chen et al, 2011). 
Genetic markers from the nuclear genome are essential 
for us to further validate this conclusion and to analyze 
other behaviors of this animal, e.g., female-biased 
dispersal and gene flow between different populations. 
Microsatellites are short tandem repeats (STR) in 
the genome and have been broadly used for evaluating 
genetic diversity, population affinity, gene flow and 
population microdifferentiation (Goldstein et al, 1999; 
Schl6etterer & Pemberton, 1994; Waits et al, 2000), as 
well as individual identification and parentage tests (Behl 
et al, 2002; Butler, 2005). So far, there have been only 
sporadic reports for characterizing microstatellites in the 
tree shrew (Munshi-South, 2008; Munshi-South & 
Wilkinson, 2006; Srikwan et al, 2002). Amplification of 
the Chinese tree shrew genomic DNA using 
microsatellites developed from other species of 
Scandentia revealed some problems, such as lack of 
polymorphism (He et al, 2009) and had non-specific 
amplification in our samples (data not shown). Therefore, 
it is necessary to develop a set of species-specific 
microsatellite markers based on the available genome of 
the tree shrew. In this study, we successfully retrieved 12 
polymorphic microsatellite markers from Chinese tree 
shrews and characterized them in a population of 117 
individuals. Our results showed that this set of 
microsatellite loci had sufficient data to measure 
population heterozygosity and individual discrimination. 


MATERIALS AND METHODS 


Experimental animals 

A total of 117 Chinese tree shrews inhabiting the 
suburb of Kunming, Yunnan Province, China, were 
collected and raised at the Experimental Animal Core 
Facility of the Kunming Institute of Zoology, Chinese 
Academy of Sciences. All sampling procedures were 
approved by the Institute’s Institutional Animal Care and 
Use Committee. Genomic DNA was extracted from the 
blood, ear and/or muscle tissues of each individual by 
using the genomic DNA extraction kit (Tiangen, Beijing). 
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Microsatellite identification and amplification 

Based on the genomic DNA sequences of the 
northern tree shrew (Tupaia belangeri) provided by the 
Ensembl database (ftp://ftp.ensembl.org/), 20 loci were 
chosen with a core repeat unit of 4—6 nucleotides. Primer 
pairs were designed using primer premier 5. PCR 
amplification condition for each locus contained a pre- 
denaturation cycle at 94 ‘C for 5 min, followed by 30 
cycles of denaturation at 94 °C for 30 s, annealing at 
55 °C for 30 s, and extension at 72 °C for 30 s, followed 
by a final extension cycle at 72 °C for 5 min. 
Amplification was performed in a total volume of 20 uL, 
which contained 10xTaq polymerase reaction buffer 
(100 mmol/L Tris-HCl, pH 8.3; 500 mmol/L KCI; 15 
mmol/L MgCl»), 0.2 mmol/L dNTPs, 10 umol/L of each 
forward and reverse primer, 0.5 U Taq DNA polymerase 
(TaKaRa, Dalian, China), and 50 ng genomic DNA. PCR 
products were electrophoresed on 8% polyacrylamide 
gels to test successful and specific amplification. For 
those loci with no amplification or with non-specific 
bands, we optimized the reaction conditions by 
increasing or decreasing the annealing temperature, or 
adjusting the concentration of Mg” to achieve the best 
efficiency and specificity. We validated the authentic 
sequencing of each amplified locus by TA-cloning 
sequencing. In brief, PCR product was ligated into PMD 
19-T vector (TaKaRa, Dalian, China), followed by 
transforming DH5a competent cells (Tiangen, Beijing) 
according to the manufacturer’s instructions. We picked 
up three to five plasmids with the right inserts for 
sequencing. 

After the pretest for primer pairs and the opti- 
mization for amplification, we labeled the forward 
primer of each locus by including 6-FAM (Carboxy- 
fluorescein) at the 5' end of the primer. PCR 
amplification for each locus was performed using the 
above optimized conditions. PCR product was diluted 
with double-distilled water to achieve a concentration of 
1 ng/uL. We mixed 1 uL diluted PCR product with 0.1 
uL Genescan’™-500 LIZ” (Applied Biosystems) and Hi- 
Di™ Formamide (Applied Biosystems), and added 
ddH,0 to make a 10 uL volume of cocktail. The mixture 
was denatured at 95 °C for 5 min, then immediately 
chilled on ice. We loaded 1 uL of cocktail on an 
automated sequencer (ABI PRISM 3100, Applied 
Biosystems). Alleles were scored using program 
Genemarker V1.5 (SoftGenetics LLC, State College, PA). 


Statistical analysis 

Allele frequencies, observed heterozygosity (Ho), 
expected heterozygosity (He) and deviation from the 
Hardy-Weinberg equilibrium (HWE) based on the 
likelihood ratio test were calculated by POPgene 1.32 
(Yeh & Boyle, 1997). Polymorphism information content 
(PIC), power of discrimination (PD), power of exclusion 
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(probability of excluding relatives of the true father from 
paternity, PE) and paternity index (PI, which is the ratio 
of the probability that a genetic marker/allele the alleged 
father passed to the child to the probability that a 
randomly selected unrelated man of similar ethnic 
background could pass the allele to the child) were 
calculated by PowerStats V1.2 (Promega, Madison, WI). 
Non-exclusion probabilities (NEP), probability of 
identity (Pid), probability of identity of siblings (Pid-sib) 
and null allele frequency (F (null)) were calculated using 
Cervus 3.0 (Kalinowski et al, 2007). 


RESULTS 


Characterization of the microsatellites 
Among the 20 selected loci, only 14 could be 


successfully amplified in the Chinese tree shrew using 
the primer pairs designed in this study (Table 1). The 
sequence authenticity of each locus was verified by 
sequencing and the microsatellite sequences were 
deposited in GenBank under accession numbers 
JQ173849-JQ173882 (Figure 1). There are some 
mutations in the repeat motif (e.g., from AAG to GAG in 
loci TB9 and TB14), besides the change of the number of 
repeats. The remaining 6 loci either failed to be amplified 
or with apparent non-specific amplification (Table S1). 
To test whether the 14 loci were polymorphic in the 
Chinese tree shrew population, we screened 117 
individuals by using 5'-fluorescently labeled primers. 
Two loci, TB2 and TB10, were non-polymorphic in the 
tested samples. However, a comparison of the sequenced 


Table 1 Information of 14 primer pairs for amplifying microsatellite loci in Chinese tree shrew 
(Tupaia belangeri chinensis) and PCR condition 


Locus Primer sequence (5'-3') 


F:ATCAGAATCTGGTTTCAAAGGT 
TB1 214-224 
R: GCACACCATGATGTAGCTGT 


F:AATATGCTGGTAGGACCAGAAAC 
TB2 185° 
R: ATTGGGTTGTTGGGAAAGTC 


F:TCTGGACCTAGCAACTTGCTTTC 
TB3 255-271 
R: TTCTTCACCAGTGTTGAGCATAA 


F:AGACAGAATGCAAGAAATCAC 
TB6 410-431 
R:ATGTGCAATGTAATAGTTCCAG 


F:CTCTGACTTCAAGACCTGGAAC 
TB8 449-473 
R:CTTGAAAGCAAAGGCAGTACAA 


F:CCTTGTGGCTTGAAAAGTATGCT 
TB9 345-389 
R:AGCGAGAATCATCCTCCAGGG 


min F:ATTGGAGACCACTCTACCCTTTA joe 
R: TGGATTCCCACATAGTCTACCG 


F:GTTTTTCCTCTCTCATGTACTC 
TB12 360-404 
R:GTCCACATCTTTGCATATCTCC 


F:GACCACGGTTCAAGTCCTAG 
TB14 415-479 
R:GATCACAGGGCTGACCTATTA 


F:TGCGGTATCATTGTGGTGGTTCA 
TB15 237-256 
R:GGGCACGTGTGCACATACCTCA 


F:TCCTTAGTATGGACTTGTTTGC 
TB16 206-256 
R:AGGATCACAGGGCAGGGAGACT 


F:ATTACGTCATCCGCTGGTCC 

TB17 290-315 
R:GCAGATTCAGACTTGCCTATCC 
F:TTGACCTAAGAGGTAAATACTG 

TB18 426-498 
R:AAGGTGAGTTTACAGCATAGAC 
F:TGATCTCCAAAGGCTTCTCCAT 


R:GAACCTTCCCGTATGGGTCAGT 


209-229 


Product (bp) 


Annealing temperature ° Mg” (mM) Repeat element 

56 °C 1.5 (AAGAAA)n 
‘ces, S870 fradaycies 15 (AAG)n(GAGAAGn 
Cela S8C forsdeytes | TAA 

56 C 1 (AAAAC)n...(TTA)n 

55 °C 1 ((A/G)TC(C/T))n 
cles S8C for26eyelee, | AVIA 

50 °C 1 (ATTT)n 
(68 —60 °C}-1'C/cycle for 9 1s (T(T/C)(T/C)C)n(TC)n 
cycles, 58 °C for 26 cycles (TTTCCC)n 
oncla S8 C helsecin |5 CAIAAG a 
G a 
cycles 580 for 30 oyes | (OMON 
ma 6 6 Eo Oa 
eales 58 © for d6 eyeles 18 TTC 
(64—60 °‘C)-1°C/cycle for 5 i (AAAAC)n 


cycles, 58 °C for 30 cycles 


“non-polymorphic microsatellite markers in Chinese tree shrews; 
10% DMSO were included in the PCR reaction; 


“PCR amplifications for all loci were performed with 35 cycles. With the exception of TB1, TB6, TB8, and TB10, the other loci were amplified with a 


touch-down protocol. For instance, TB2 was amplified using a decreasing annealing temperature by 1 °C every cycle from 64 ‘C to 60 © (total 5 


cycles), then the PCR amplification was performed with an annealing temperature of 58 °C for 30 cycles. 
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Figure 1 GeneScan maps and sequences of 12 polymorphic microsatellites in Chinese tree shrew population 


fragments of these two loci with the reference sequence 
from the Ensembl database showed length 
polymorphisms (Figure S1), suggesting that these two 
loci might present polymorphisms among tree shrews 
from different regions, but not in tree shrews collected 
from Kunming. This result was in accordance with our 
previous studies of mtDNA sequence variation in the 
Chinese tree shrew (Chen et al, 2011; Xu et al, 2012), in 
which we observed some nucleotide differences between 
the Chinese tree shrew from Kunming and the reported 
northern tree shrew outside China, though both were 
grouped as Tupaia belangeri. 


Chinese tree shrew population has a high microsatel- 
lite polymorphism 

Tables 2 and 3 list the estimated parameters for each 
of the 12 polymorphic microsatellites in Chinese tree 
shrews. The number of alleles per locus ranged from 4 to 
20, with an average of 8.75 alleles. The effective alleles 
ranged from 1.19 (TBI locus) to 10.47 (TB18 locus), and 
had an average of 4.05. The mean values of observed 
heterozygosity and expected heterozygosity were 0.566 
and 0.616, respectively. Loci TB1 and TB15 had a 
relatively low level of polymorphism (PIC value<0.25), 
TB3 and TB6 were moderately polymorphic (0.25<PIC 
value<0.5), and the other loci were highly polymorphic 
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(PIC value>0.5) (Table 2). The high PIC value and 
observed heterozygosity indicated a high genetic 
diversity of the analyzed tree shrew population, which 


Table 2 Genetic diversity measures for 12 polymorphic 
microsatellite loci tested in 117 Chinese tree shrews 


Locus Na Ne Ho He PIC 


TBI 4 1.193 0.103 0.107 0.12 
TB3 6 1.786 0.436 0.442 0.40 
TB6 6 2.226 0.537 0.553 0.46 
TB8 7 3.733 0.675 0.735 0.69 
TB9 7 3.256 0.709 0.696 0.65 
TB12 20 7.517 0.607 0.871 0.86 
TB14 17 8.064 0.829 0.880 0.86 
TB15 4 1.213 0.137 0.177 0.17 
TB16 7 3.041 0.675 0.674 0.62 
TB17 5 3.247 0.641 0.695 0.64 
TB18 18 10.47 0.838 0.908 0.90 
TB20 4 2.897 0.607 0.658 0.59 
Average 8.75 4.054 0.566 0.616 0.58 


Na: observed number of alleles; Ne: effective number of alleles; Ho: 
observed heterozygosity; He: expected heterozygosity; PIC: 


polymorphism information content. 
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was consistent with the previous mtDNA analysis (Chen 
et al, 2011). A likelihood ratio test for concordance with 
HWE revealed that only loci TB12, TB17 and TB20 
deviated from HWE (P<0.05), but none was significant 
after applying a Bonferroni correction. 

To further evaluate the potential use of the 12 
microsatellite markers, we calculated some forensic 
efficiency parameters based on the tested samples. The 
linkage disequilibrium was tested by the online SHEsis 
program (Shi & He, 2005), but no significant linkage 
between loci was detected (Table S2). The power of 
discrimination (PD) was from 0.218 (TB1) to 0.971 
(TB18), and the cumulative PD in this study was as high 
as 99.99999998%. Waits et al (2001) recommend a 
reasonable probability of identifying individuals in a 
given region using dominant and codominant markers as 
0.01—0.0001. In this study, the values of combined Pid 
and combined Pid-sib were 4.89x10"' and 0.00001, 
respectively. These values are far smaller than 0.01, 
which suggests that these 12 loci had sufficient data for 
individual identification. 

As testing for more markers would increase the cost, 


time and the risk of analytical errors (Waits & Paetkau, 
2005), it was necessary to choose the minimum number 
of markers that fulfill the aim of individual identification. 
We eliminated the four less polymorphic loci (TB6, TB3, 
TB15 and TB1) but still obtained considerably high 
values for combined PD (CPD=99.99999955%), 
combined Pid (CPid=8.61x10"°), and combined Pid-sib 
(CPid-sib=5.28* 10°). It is thus evident that the remaining 
8 loci are sufficient for individual discrimination. 

For parentage testing, parameters such as average 
non-exclusion probability (NE-1P), average non- 
exclusion probability for one candidate parent given the 
genotype of a known parent of the opposite sex (NE-2P), 
average non-exclusion probability for a candidate parent 
pair (NE-PP), power of exclusion (PE) and paternity 
index (PI) were widely used in the field. The Chinese 
tree shrew population had a NE-1P value of 0.322—0.994, 
a NE-2P value of 0.192—0.945, a NE-PP value of 
0.059—0.897, a PE value of 0.668—0.012, and a PI value 
of 3.05—0.57, respectively (Table 3). Excluding the four 
loci (TB3, TBI5, TBl, and TB6) with few 
polymorphisms, we still obtained sufficient data for 
parentage testing. 


Table 3 Individual identification and parentage testing parameters for 12 polymorphic microsatellite markers in a 


population of Chinese tree shrew (n=117) 





Locus F (null) PI 





TB18 0.034 3.05 3.05 0.322 0.32 0.192 0.19 0.059 


CPI NE-1P CNEIP NE-2P CNE2P NE-PP CNEPP Pid CPid 


Pid-sib CPid-sib PD CPD PE CPE 


0.017 0.02 0.302 0.30 0.971 0.971 0.668 0.668 


TB14 0.027 2.93 8.94 0.401 0.13 0.249 4.78x10° 0.095 5.54x10° 0.028 4.64x107* 0.319 9.63107 0.966 0.9990140000 0.654 0.147 


TB9 -0.006 1.72 15.37 0.72 9.29x107 0.544 2.60x107 0.358 1.98x10° 0.138 6.43x10° 0.438 4.22x107 0.849 0.9998511140 0.443 0.337 
TB8 0.04 1.58 24.29 0.673 6.25107 0.496 1.29x10° 0.309 6.14x10% 0.113 7.23x10° 0.412 1.74107 0.886 0.9999830270 0.404 0.605 


TB16 —0.002 1.54 37.40 0.74 4.62x10° 0.57 7.35x10° 0.386 2.37x10% 0.157 1.14x10° 0.454 7.89x10° 0.837 0.9999972334 0.391 0.78 


TB17 0.044 1.39 51.99 0.727 3.36107 0.559 4.11x10° 0.38 9.00x10° 0.148 1.68x10°" 0.441 3.48x10° 0.864 0.9999996237 0.343 0.846 


TB12 0.19 1.27 66.02 0.404 1.35x107 0.252 1.04x10° 0.089 8.01x10° 0.028 4.67x10° 0.032 1.13107 0.935 0.9999999755 0.299 0.947 
TB20 0.046 1.27 83.85 0.776 1.05x10° 0.621 6.43x10% 0.461 3.69x10° 0.184 8.61x10'° 0.469 5.28x10° 0.816 0.9999999955 0.299 0.948 


TB6 0.011 1.08 90.56 0.845 8.89x10° 0.734 4.72x10% 0.599 2.21x10° 0.292 2.52x10'° 0.548 2.89x10° 0.71 0.9999999987 0.223 0.968 


TB3 —0.009 0.89 80.60 0.899 7.99x10° 0.763 3.60x10% 0.616 1.36x10° 0.354 8.90x10'! 0.618 1.79x10° 0.653 0.9999999996 0.137 0.979 
TB15 0.115 0.58 46.75 0.984 7.86x10° 0.908 3.27x107 0.83 1.13x10° 0.685 6.10x10'' 0.833 1.49x10° 0.291 0.9999999997 0.015 0.993 


TB1 0.048 0.57 26.64 0.994 7.82x10° 0.945 3.09x10% 0.897 1.01x10° 0.801 4.89x10'! 0.897 1.34x10° 0.218 0.9999999998 0.012 0.995 


F (null): null allele frequency; NE-1P: Average non-exclusion probability for one candidate parent; CNE1P: combined NE-1P; NE-2P: Average non- 
exclusion probability for one candidate parent given the genotype of a known parent of the opposite sex; CNE2P: combined NE-2P; NE-PP: Average 


non-exclusion probability for a candidate parent pair; CNEPP: combined NE-PP; Pid: probability of identity, the probability of that two individual 


share the same genotypes using all loci; CPid: combined Pid; Pid-sib: the probability of identity of siblings; CPid-sib: combined Pid-sib; PD: power 


of discrimination; PE: power of exclusion; CPD: combined PD; PI: paternity index; CPI: combined PI. 


The null allele frequency was also calculated to 
evaluate the efficiency of this set of microsatellite markers 
because of its influence on population analysis (Chapuis 
& Estoup, 2007). Null allele is any microsatellite allele at 
a microsatellite locus that consistently fails to amplify to 
detected levels in the PCR assays, and it has no significant 
effect on parentage analysis when the frequency is less 
than 0.2 (Dakin & Avise, 2004). In our study, all the null 
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allele frequencies were less than 0.2 (Table 3). Therefore, 
these markers can be used for population analysis of tree 
shrews, at least in parentage analysis with little influence 
on the average exclusion probability. 


DISCCUSSION 


Tree shrews have received wide attention in recent 
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years for their potential use in biomedical research and 
drug discovery (Cao et al, 2003; Fuchs & Corbach-Sohle, 
2010; Peng et al, 1991). Despite decades of cultivation 
efforts, no inbred strain has been established for the tree 
shrew, which hindered a wide use of this animal in 
biomedical research; this might account for the 
inconsistency of previous results concerning animal 
models. Thus establishing some inbred lines of tree 
shrews is crucial, and we have launched an ambitious 
project to fulfill this goal. 

To facilitate our current inbreeding program for the 
Chinese tree shrew and to evaluate the genetic diversity 
of wild tree shrews, we needed a set of microsatellite 
markers that had sufficient data for individual 
discrimination and parentage testing, as well as for the 
quantification of population heterogeneity. In this study, 
by referring to the available tree shrew genome sequence 
(albeit the coverage was only 2x), we were able to obtain 
12 polymorphic loci in the Chinese tree shrew. 
Evaluation of different parameters related to genetic 
diversity, individual identification and parentage testing 
showed that eight of these loci had sufficient data for 
lineage tracing and measurement of diversity. 

Compared with these reported microsatellite 
markers, which were developed for T. glias and T. minor 
(Munshi-South & Wilkinson, 2006; Srikwan et al, 2002), 
our newly designed set based on the genome sequence of 
T. belangeri was specifically optimized for Chinese tree 
shrews. As we had no samples of other species in Tupaia, 
we did not know whether our primer pairs would work 
for these species. As Chinese tree shrews were proposed 
to contain six subspecies (Wang, 1987), it would be 
rewarding to screen the genetic diversity of these 
different subspecies, both for conservation of genetic 
resources and for selecting a proper founder to establish 
the inbred lines. 


References 


Arnason U, Adegoke JA, Bodin K, Born EW, Esa YB, Gullberg A, 
Nilsson M, Short RV, Xu X, Janke A. 2002. Mammalian mitogenomic 
relationships and the root of the eutherian tree. Proc Natl Acad Sci USA, 
99(12): 8151-8156. 


Behl R, Sheoran N, Behl J, Tantia MS, Vijh RK. 2002. Microsatellite 
sequences of mammals and their applications in genome analysis in 
pigs: A review. Asian Austral J Anim, 15(12): 1822-1830. 


Butler JM. 2005. Forensic DNA Typing: Biology, Technology and 
Genetics of STR Markers. New York: Elsevier Academic Press. 


Cao J, Yang EB, Su JJ, Li Y, Chow P. 2003. The tree shrews: adjuncts 
and alternatives to primates as models for biomedical research. J Med 
Primatol, 32(3): 123-130. 


Chapuis MP, Estoup A. 2007. Microsatellite null alleles and estimation 
of population differentiation. Mol Biol Evol, 24(3): 621-631. 


Chen SY, Xu L, Lu LB, Yao YG. 2011. Genetic diversity and 


Kunming Institute of Zoology (CAS), China Zoological Society 


Microsatellite markers have long been used in 
genetic analyses for both wild populations and/or 
domestic animals. In general, wild populations have a 
higher genetic diversity compared to the inbred strain or 
domestic breeds. Analysis of the inbred C57BL/6J mouse 
strain using microsatellite markers revealed a lower 
expected heterozygosity (<0.5) (Niu & Liang, 2009), 
but the expected heterozygosity of the wild house mouse 
(Mus musculus) reached 0.62 in population from 
Cameroon and 0.82 in population from Germany 
(Thomas et al, 2005). Based on the analyses of 19 
microsatellite loci, three pure breeds of dog (Greyhounds, 
Labradors and German Shepherds) showed expected low 
heterozygosities of 0.357, 0.481 and 0.431, respectively 
(Zajc et al, 1997). In this study, the expected 
heterozygosity of the Chinese tree shrew was 0.616, 
lower than that of the plain treeshrew (T. longipes) (0.74) 
but higher than that of the large treeshrew (T. tana) (0.58) 
distributed in Sabah, Malaysia (Munshi-South & 
Wilkinson, 2006). As our population size was considerably 
large (n=117), we speculated that the observed 
heterozygosity might reflect the true situation of wild 
tree shrews distributed in Kunming, Yunnan Province. 

In short, we characterized a set of 12 polymorphic 
microsatellite markers for the Chinese tree shrew, which 
were confirmed to be sufficient for assessing population 
genetic structure and individual discrimination of this 
animal. An analysis of 117 wild individuals showed that 
the Chinese tree shrew had a considerably high 
heterozygosity. We hope that these markers will provide 
essential help for us to advance the inbreeding project of 
the Chinese tree shrew. 
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Table S1 Information of the 6 pairs of primers that could not be amplified or had 
non-specific amplification 


Locus 


TB4 


TB5 


TB7 


TB11 


TB13 


Primer sequence (5'-3') 


F:GTCCTAGACCTGGCAGCTTGCT 
R:TTGCTATGTGACCATCCCTGTG 
F:GGGTAGATACCCAGGAGTAGGAT 
R:AGAGATTGCAGTGTGATCATGC 
F: TAAGGACCTGAAAGATGTTTGC 
R:TCCCATGACTGACTGTGGTTTG 
F:TTTGTACTTTATGTAGTTACTC 
R:TTGAAAACACATCAGGCAC 
F:TGGCTTAGACTCCTGCTCGGTG 
R:TTACCACTGCACCATCCTGACG 
F:GTAGCTTGAAAATCATGACAGGC 


R:GAATGCTAGAAACACTCAGTGGC 


Primer 


length 
(bp) 


GC% 


Predicted 
product 
(bp) 

400 

426 

582 


450 


367 


Repeat element 


(AGAT)n 


(CCTTT)n (CCTTC)n 


(TTTC)n (TTTCTC)n 


(TTTC)n(TTCC)n 


(AAAG)n(GAAG)n 


(GAAA)n 


Liu and Yao 2013 Zool Res 34(E2): E62-E68 


Table S2 Linkage disequilibrium of 12 polymorphic microsatellite markers in Chinese tree shrew populations 


TB1 
TB3 
TB6 
TB8 
TB9 
TB12 
TB14 
TB15 
TB16 
TB17 
TB18 
TB20 


TBI 


KK 


0.026 
0.001 
0.012 
0.007 
0.005 
0.006 
0.001 
0.018 
0.022 
0.007 
0.001 


TB3 
0.389 


KK 
0.017 
0.018 
0.017 
0.006 
0.007 
0.01 
0.02 
0.01 
0.013 
0.008 


TB6 
0.145 
0.217 
KK 
0.007 
0.014 
0.01 
0.01 
0.002 
0.013 
0.009 
0.01 
0.004 


TB8 
0.443 
0.304 
0.17 
KK 
0.016 
0.007 
0.014 
0.01 
0.012 
0.017 
0.007 
0.017 


TB9 
0.346 
0.308 
0.248 
0.297 
KK 
0.007 
0.008 
0.023 
0.019 
0.008 
0.014 
0.005 


TB12 
0.493 
0.336 
0.393 
0.352 
0.346 
KK 
0.01 
0.013 
0.009 
0.006 
0.008 
0.008 


TB14 
0.481 
0.289 
0.304 
0.357 
0.3 
0.456 
KK 
0.009 
0.008 
0.009 
0.009 
0.011 


TB15 
0.332 
0.207 
0.134 
0.352 
0.637 
0.519 
0.459 
KK 

0.021 
0.002 
0.008 
0.006 


TB16 
0.453 
0.254 
0.235 
0.264 
0.302 
0.358 
0.302 
0.484 
KK 

0.003 
0.011 
0.011 


TB17 
0.461 
0.197 
0.176 
0.268 
0.205 
0.306 
0.302 
0.157 
0.122 
KK 

0.012 
0.014 


TB18 
0.572 
0.387 
0.346 
0.349 
0.408 
0.491 
0.47 

0.522 
0.358 
0.368 


KK 


0.012 


TB20 
0.142 
0.17 
0.126 
0.243 
0.119 
0.3 
0.283 
0.236 
0.2 
0.229 
0.349 


KK 


The upper and lower triangular matrixes present the D' values and the 7° values, respectively. The values were estimated based on 117 Chinese 


tree shrews. 
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(A) TB2R AATATGCTGGTAGGACCAGAAACAAAGGAAGGAACTAAATAAAGATTITTTTAAAAGAAG 
TB2C AATATGCTGGTAGGACCAGAAACAAAGGAAGGAACTAAATAAAGATTTTTTTAAAAGAAG 
TB2R AAGAAGAAGAAGAAGAAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAG 
TB2C AAGAAGA——---—-----------—-------—--—------—------------ AG 
TB2R AAGGAGAAGGAGAAGGAGAAGAAGGAGAAGAAGAGAATGAACATATTAGGTCATCAGGAG 
TB2C AAGGAGAAGGAGAAGGAGAAGAAGGAGAAGAAGAGAATGAACATATTAGGTCATCAGGAG 
TB2R ATTCATGACATATTTTCCTATGCAAAAATGCATCATGACTTTCCCAACAACCCAAT 
TB2C ATTCATGACATATTTTCCTATGCAAAAATGCATCATGACTTTCCCAACAACCCAAT 

(B) Bior GGACTTTATGCCTCATATCCAACCCTCTTGGGTATTTATTTATTTATTTATTTATTTATT 
TB10C GGACTTTATGCCTCATATCCAACCCTCTIGGGTCTTTATTTATTTATTTATTTATTTATT 
TBIOR TATTTATTTATITATTTATTITATITATTTTTGGAGAGGTTTTCGCTACAGGGCGCATGCG 
TB10C T---------------—----------- TTGGAGAGGTTTTCGCTACAGGGCGCGTGCG 
TB10R CTCGCCAGACAGAGAGACATAGAGACAGAAAGCGAGAAATCACCCCCCAGAGGGGATGGT 
TB10C CTCGCCAGACAGAGAGACATAGAGACAGAAAGCGAGAAATCACCCCCCAGAGGGGATGGT 
TB10R CACGCCGCAATTCCTTCTCCTCCATTCTCCTCTCGGGATCTTGGGGACGCACGCACGCGC 
TB10C CACGCCGCAATTCCTTCTCCTCCATTCTCCTCTCGGGATCTTGGGGACGCACGCACGCGC 


Figure S1. Alignment of sequences of loci TB2 (A) and TB10 (B) in Chinese tree 
shrew with the tree shrew reference from the Ensembl database. There are differences 
for the number of repeat motif in Chinese tree shrews (locus TB2: TB2C; TB 10: 
TB10C) and the tree shrew from the Ensembl database. 


